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Abstract 



On the scales of galaxies and beyond there is evidence for unseen dark matter. 
In this paper we find the experimental limits to the density of dark matter 
bound in the solar system by studying its effect upon planetary motion. 
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1. INTRODUCTION 



According to Newton's inverse square force law, the circular speed around an isolated 
object of mass M should be 

[MG 

In disk galaxies we do, however, observe that the circular speeds are approximately indepen- 
dent of r at large distances. The standard explanation is that this is due to halos of unseen 
matter that makes up around 90% of the total mass of the galaxies (Tremaine 1992). The 
same pattern repeats itself on larger and larger scales, until we reach the cosmic scales where 
a baryonic density compatible with successful big bang nucleosynthesis is less than 10% of 
the density predicted by inflation, i.e. the critical density. 

The flat rotation curves of galaxies, taken at face value, imply that the effective gravi- 
tational force follows a 1/r law at large scales. This could either be due to dark matter or 
to a departure from Newtonian dynamics at small accelerations (Milgrom 1983; Bekenstein 
1992) or large scales (Sanders 1990). An effective gravitational acceleration law of the form 

9 = 

r 

at small accelerations a < do has been reported (Kent 1987; Milgrom 1988; Begeman, 
Broeils, & Sanders 1991) to be successful in reproducing the observations of galactic sys- 
tems.[]The constant a has been determined by studies of galaxy rotation curves and its 
value has been found to be a « 10~ 8 cm s~ 2 . As noted by Milgrom (1983), this value of 
a « cH . 

With such a 1/r force law the circular speed would approach v c = (GMao) 1 ^. If the 
luminosity L of a galaxy is proportional to its mass M, then this relation would explain the 
infrared Tully-Fisher law (Tully & Fisher 1977) which states that circular speeds in galaxies 
scale as v c oc 

The theoretical underpinning for the 1/r effective force law is not yet firmly established. 
It might be due to a modification of gravity along the lines of Milgrom (1983), but it seems 
to be difficult to construct a viable relativistic theory of this kind (Zhytnikov and Nester, 
1994). Accordingly, the standard view is that the effective galactic 1/r force law is caused by 
dark matter. At this point, it is worth mentioning that a large-distance force law of this type 
can be reproduced within standard general relativity theory with a very simple, but perhaps 
unrealistic, matter source (Soleng 1993, 1995). Our key point is that general relativity is 
quite capable of explaining the observed gravitational properties of the universe provided 
we give it the right input. Most likely the dark matter is a mixture of several components, 
such as weakly interacting particles, black holes, brown dwarfs, neutron stars, as well as 
energy stored in high-frequency oscillation of Newton's gravitational coupling (Accetta & 
Steinhardt 1991; Steinhardt & Will 1994). Whatever the origin of the 1/r force law might 



However, not without debate (Lake 1989; Milgrom 1991). 



2 



be, its reported experimental success forces us to take it seriously. Accordingly, we think 
that it is particularly important to compare the densities of dark matter inferred from large 
scale dynamics with experimental limits from local tests. If dark matter exists in the form 
of microscopic objects, one would expect that this unknown form of energy penetrates into 
galaxies and also enter the solar system. 

Braginsky, Gurevich, and Zybin (1992) have studied the effect of dark matter bound in 
the galaxy but unbound to the solar system. Such unbound dark matter would produce an 
anisotropy in the gravitational background of the solar system. The resulting tidal forces 
induce an additional perihelion precession. Assuming = 0.3 GeV/cm 3 = 5.4 x 10~ 25 g/cm 3 
they computed the magnitude of the resulting secular orbit distortion. The effect may be 
observed by reasonable improvements of present observational techniques (Klioner and Soffel 
1993; Braginsky 1994). A possible influence of dark matter on the Earth-Moon system has 
been considered by Nordtvedt (1994) and by Nordtvedt, Miiller, and Soffel (1995). 

In this paper we focus on a dark matter model in which the density of dark matter 
varies so slowly within the solar system that it can considered constant. This is a reasonable 
assumption if dark matter in the solar system really is in the much deeper potential of the 
galaxy with the Sun causing only a local density perturbation in the galactic dark matter 
background. It will also be assumed that the equation of state of the unseen matter is 
almost dust-like, that is, the pressure will be assumed to be much less than the energy- 
density. Based on this model we calculate an upper limit to the density of dark matter by 
considering its effect upon the perihelion precession of the planets. We have also carried out 
similar computations with a dark matter density proportional to 1/r 4 and 1/r 2 using the 
results in Soleng (1994) and (1995), respectively. The corresponding experimental bounds 
do not vary more than one order of magnitude. This weak dependence on the distribution 
function corroborates the claim of Anderson et al. (1989) and should be expected because 
(to lowest order) the perihelion precession caused by dark matter is given by the integrated 
dark matter mass at a given orbital radius. 

2. SOLAR SYSTEM WITH DARK MATTER 

In order to study the gravitational effects of hypothetical dark matter on planetary 
motion, we need a solution of Einstein's field equations for a static, spherically symmetric 
space-time and a given distribution of dark matter. The line-element for a static, spherically 
symmetric gravitational field can in general be written asQ 

dg 2 = _ e 2Kr) dt 2 + e 2X(r) dr 2 + ^^2 (1) 

We shall assume that the dark matter has a constant density po (within the solar system). 
At a surface where the dark matter pressure equals the galactic dark matter pressure pc ~ 
10~ 7 pG (characterized by velocities of 220km/s), we match the gravitational field of the solar 
system to the exterior field of the galaxy. We shall assume that g a ~ —1 at this distance 



2 We employ geometrized units with G = c = 1. 
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(this assumption is always used in local gravitational problems). Then the tt-component of 
the field equations takes the form 



d 
dr 



r 1 



-2A 



inp r 



Integrating with a spherical mass M (the Sun) at r = gives 

2M 8tt „ 



-2A 



1 



-Por 



r 3 

The Tolman-Oppenheimer-Volkov (TOV) equation takes the form 

dp _ M + ^-p r 3 + inpr 3 

~dr ~ ~ (P ° + P K 2 - 2r(M + f p r 3 ) ' 



(2) 



where p is the pressure of the dark matter. Using that the mass of dark matter in the 
solar system is much less than the mass of the Sun, we can neglect the last term in the 
denominator of equation (Q). Also we assume that the dark matter is non-relativistic, i.e. 
that p <C po, and accordingly we can neglect the last term in the numerator. Then the TOV 
equation reduces to 



dp , , ,M + fp r 3 
d^ ~ ~ {P0 + P} r(r - 2M) ' 



Integration leads to 



P 



1 _ 2M V2M 



-16A/ 2 7rp /3 



exp 



~p Q (4Mr + r 2 ) 



Po, 



where K\ is a constant of integration. With the assumptions that M 2 po <S 1 and M < r, 
this equation simplifies to 



P 



-e s 



-Pa- 



(3a) 



The constant K\ can be determined as follows: let po be the pressure at r = r sur f where r sur f 
is the surface of the central mass (the Sun in our case). Within our approximations, 



Ki = Po+Po- 



(3b) 



The relativistic equation of hydrostatic equilibrium, T M ^ ;i , = 0, where T^ v is the energy- 
momentum tensor of the dark matter, leads to 



dp , dp: 
dr dr 



Integration yields 



9tt = -e 



2/i 



Po+Po 
,Po+P, 



where K 2 is a new integration constant. Inserting the pressure from equations ([3D leads to 



2M 



gtt = ~K 2 \ l - —J e 3 

The constant K 2 can be determined by demanding that g u = — 1 at the surface r = r matc h 
where the pressure equals the galactic dark matter pressure Pq. This radius is given by 



'"match 

and hence 



_3_ ln Po+Po 
2vrp \Po + Pg, 



1/2 



(4) 



r - i ^ 2 \ [P0+ PG 

K 2 = exp -— p r match - 



2 

o <"' uiatcii / i ) ~ 1 

3 / V Po + Po J 

according to the assumption that p <C p . 

Our model of the dark matter filled space—time in the solar system is thus represented 
by the approximate line-element 

, 9 / 2M 47r 9 \ , o dr 2 9 ~ 

<fc 2 = - (l - — + - Po r 2 ) dt 2 + 8 _^ por2 + rW (5) 

valid outside the Sun and for r <C (M/po) 1//3 - It should be noted that we only require the 
validity of this expression inside the solar system, where ^fp$r 2 is a small number. At larger 
scales where the dark matter pressure equals its galactic value, the gravitational field is of 
course determined by the mass distribution of the galaxy. 



3. PERIHELION PRECESSION 

The Lagrange function for a test particle moving in the 9 = 71/ 2 plane in the geometry 
specified by equation (g), is 

21 = - 1 1 - — + T"» r ) ' + 1 - S - &p + r 41 

where a dot denotes differentiation with respect to the proper time r of the particle (a 
planet). The momenta p M = dL/dx^ 1 are 

/ 2M 4?r 2 \ . 
Pt—\1 h ~~^~Po r J t = constant (6a) 

/ 2M 8tt ,\ _1 
Pr = (1 —por J r (6b) 

= r 2 cp = constant. (6c) 
Using the normalization of the momenta g liv p tl p u = — 1, we get 

9tt 9rr 
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To first order in M and po, one finds 

. 2 / 2M 8tt 2 \( p 2 A 2 2 2 

r = \1 - — - — Po r J I -1 - I + p t - 47rp tPo r . (7) 

Inserting the new radial variable u = 1/r and noting from equation (|6cD that f = 
—p^(du/d(j)), brings equation ([/]) into the form of an orbit-equation. Differentiation of the 
resulting expression, and using that p t « 1 for planets in the solar system, leads to 



" - -j- + 3Mu + — — 

pj 3 u 



+ M + 3M M 2 + --gL (8) 



The planetary orbits are nearly circular, and we can treat the perihelion precession as a 
perturbation from the circular solution, u = uq, where 

M 9 47T po 

«o = -^ + 3M«S + — -gp (9) 

Substituting uq(1 + e) for u into equation (H) with e < 1, and using equation (§) and 
calculating to first order in e, we get 

= 



The Einstein precession coming from the solar mass M = M & is 

A0 O = 6vrM M O - (10) 
In addition, there is a dark matter induced precession 



dark = -47T 2 -^2". (11) 



Equations fllUl) and QTTD imply 



= Tm^4 a *- (12) 

Let the observed non- Newtonian perihelion precession be denoted by A0 o b s , and its un- 
certainty by 50obs- Since A0 o b s = A0 O within the uncertainty, the dark matter perihelion 
precession and the Einstein term are related by |A0dark| < |#0obs|- This, together with 
equation fli~2|), gives 



3 %, „2„.5l^ 

To lowest order p 2 , = M/u . Thus 



<- ~ n/ r 2 5 0< Pobs 



< 3 GMq \5<p ohs \ 
P ° " 2tt c 2 4 A0 O 
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where we have inserted Newton's constant G and the speed of light c in order to simplify 
the numerical calculations. The perihelion precession of the asteroid Icarus is known with 
about 8 % accuracy. Its distance from the Sun is 1.076 A.U. With M & = 2 x 10 33 g and 
r = 1.076 A.U.= 1.61 x 10 13 cm we obtain 

p £ 1.8 x 10" 16 g/cm 3 . (13) 

This value is about seven orders of magnitude above the mean galactic mass density, and it 
shows that measurements of the perihelion precession of the planets do not put strict limits 
on the density of bound dark matter. Note that with this density, the mass of spherically 
symmetric bound dark matter of constant density ( [IB] ) within the orbit of Pluto is less than 
2 x 10~ 5 M & and thus in agreement with our assumptions. Within Uranus we find a limit 
of 

M dark (ru) < 2 x 10- 6 M « 0.6 M e 

where M ffi is the Earth's mass. This bound is of the same size as the bound found by 
Anderson et al. (1989) by numerically analyzing how the orbit of Uranus would be affected 
by dark matter. Recently, the bound obtained by such methods has been strengthened to 
around 0.2M® (Anderseon et al. 1995), but a similar improvement should also be possible 
using improved perihelion precession observations. 

Let us finally check the assumption that the dark matter has a pressure which is much 
less than its energy density. The maximum pressure is at the center of the solar system. 
Now, if we require that the model is valid out to the Oort cloud at several thousand A.U. 
and that macthing to the galactic dark matter distrubution takes place here, we find that 
the pressure gradient necessary for hydrostatic equilibrium is very small. Using equation (|j) 
and matching at, say, r matc h = 5000 A.U. gives 

Po ~ PG ~ 0.0015. (14) 
Pg 

Since Pg ~ 10~ 7 po, this result also confirms that p <C po- 



4. DISCUSSION 

Dark matter in the solar system would increase the perihelion precession of the planets. 
Assuming a very simple model for dark matter having a constant density and being kept 
in hydrostatic equilibrium by a small but positive pressure, we have used the perihelion 
precession of Icarus to put a limit of p £ 1.8 x 10~ 16 g/cm 3 on the density of such dark 
matter in the solar system. This bound is much larger than the mean galacitic density, but 
we note that the bound is only weakly dependent on the dark matter distribution function 
and that some dark matter models may predict a local concentration of dark matter energy. 

If dark matter exists in the form of ordinary matter, one would expect that orbital 
perturbations, especially by Jupiter, effectively would clean the inner parts of the solar 
system of any weakly interacting dark matter (Anderson et al. 1989). In this case, our bound 
which was derived from an asteroid would be irrelevant. However, cosmological constraints 
from inflationary models and cosmic nucleosynthesis, strongly suggest that a large part of 
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dark matter is non-baryonic. We do not know what the non-baryonic dark matter is. One 
possibility (Accetta & Steinhardt 1991; Steinhardt & Will 1994) is that the dark matter 
in part is oscillation energy caused by rapid oscillations of Newton's constant. An effective 
Brans-Dicke field is a consequence of many unification schemes and also an ingredient of 
extended inflationary models (La & Steinhardt 1989). Extended inflation would drive the 
scalar field away from the minimum of its potential, and the field would then start to oscillate 
when inflation ends. This possibility only illustrates that dark matter may behave rather 
differently from normal matter and that no stone should be left unturned in the search for 
the mass which seems to make up most of the universe we inhabit. 
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